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Abstract 



The six-vertex model on an iV x iV square lattice with domain wall boundary 
conditions is considered. A Fredholm determinant representation for the partition 
function of the model is given. The kernel of the corresponding integral operator 
is of the so-called integrable type, and involves classical orthogonal polynomials. 
From this representation, a "reconstruction" formula is proposed, which expresses 
the partition function as the trace of a suitably chosen quantum operator, in the 
spirit of corner transfer matrix and vertex operator approaches to integrable spin 
models. 



1 Introduction 



The six-vertex model on a square lattice with domain wall boundary conditions 
(DWBC) was introduced in [1] and subsequently solved in [2], where a determi- 
nant formula for the partition function was obtained and proven (see also [3]). 
This model, in its inhomogeneous formulation (i.e. with the vertex weights given 
as suitable functions of the position of the vertex on the lattice) , naturally arises 
when investigating correlation functions of quantum integrable models in the 
framework of quantum inverse scattering method. In its homogeneous version, 
the model admits usual interpretation as a model of statistical mechanics with 
fixed boundary conditions, and may be seen as a variation of the original six- 
vertex model with periodic boundary conditions [4-7] , which has been for decades 
a paradigmatic one in statistical mechanics [8]. 

The DWBC version of the six-vertex model has however several non trivial 
peculiarities which render it worthy of further investigations. Firstly, the model 
with this very specific fixed boundary conditions enjoys interesting connections 
with some important issues of enumerative combinatorics, such as Alternating 
Sign Matrices [9-11] and domino tilings [12,13]. Moreover, it appears that, even 
in the thermodynamic limit, bulk quantities, such as the bulk free energy, are 
indeed sensitive to the choice of boundary conditions [14,15] (see also [16]). Fi- 
nally, the determinant representation given in [2, 3] for the partition function, 
and analogous ones recently presented for the boundary one-point correlation 
functions (polarizations) [17], are rather implicit and turn out to be too intricate 
for any further, more explicit, answer, except in very particular cases. Alterna- 
tive equivalent representations for the partition function and polarizations would 
be therefore highly desirable to address several problems, such as further Alter- 
nating Sign Matrices weighted enumerations, or extensions of the Arctic Circle 
Theorem [13] beyond the Free-Fermion point. 

The model is formulated on a square lattice with arrows on edges. The only 
admitted configurations are such that there are always two arrows pointing away 
from, and two arrows pointing into, each lattice vertex; each vertex can therefore 
be in one out of six different possible states, a Boltzmann weight Wj being assigned 
to each vertex, according to its state % (i = 1,...,6). We shall consider here 
the homogeneous version of the model, where the Boltzmann weights are site 
independent. The DWBC are imposed on the N x N square lattice by fixing the 
direction of all arrows on the boundaries as follow: the vertical arrows on the 
top and bottom of the lattice point inward, while the horizontal arrows on the 
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Figure 1: The six allowed types of vertices, their weights and one of the possible 
configurations in the model with the domain wall boundary conditions for N = 5. 



left and right sides point outward. The correspondence between the Boltzmann 
vertex weights and the arrow configurations, and a typical configuration of the 
model with DWBC, are shown in Fig. 1. 

The partition function is obtained by summing over all possible arrow config- 
urations, compatible with the imposed DWBC, each configuration being assigned 
its Boltzmann weight, given simply as the product of all the corresponding vertex 
weights: 



DWBC configurations i=l 



Here rii denotes the number of vertices in state i, i.e. with Boltzmann weight Wj, 
in each configuration, and Y^h=i n-i = N 2 . The six-vertex model with DWBC is 
usually considered with its weights invariant under inversion of all arrows, and 
thus with only three distinct weight functions, denoted as a, b and c, 



wx = w 2 = a , w 3 = w 4 = b , w 5 = w 6 = c . 
We shall use the following parametrization for the weight functions 
a = sin(A + rj) , b = sin(A — 77) , c = sin2r] . 



(1-2) 



(1.3) 



In terms of this parametrization the result of Ref. [3] for the partition function 
reads 



Z N = [sin(A -; ) i Sin(A + ^ )]jV2 det N H 

n (k\y 

k=l 



(1-4) 
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where H is an N x N Hankel matrix, with entries 

„ d j+k sin 2r] 

jk d\i+ k sin(A - 77) sin(A + 77) ' 1 ' 

Here and in the following we use the convention that indices of N x N matrices 
run over the values j, k — 0, 1, . . . , N— 1. Formula (jl.4|) for the partition function, 
which was originally obtained within the Quantum Inverse Scattering Method, 
will be referred to as the Hankel determinant representation. 

The purpose of the present paper is to give some other equivalent representa- 
tions for the partition function. The emphasis will be made on the representation 
in terms of Fredholm determinant of some linear integral operator of integrable 
type, in the sense of paper [18]. Representations of such type have been proven to 
be powerful tools in many areas of mathematical physics, ranging from the the- 
ory of Random Matrices to the asymptotics of Orthogonal Polynomials. Among 
them is the theory of correlation functions of quantum integrable models [19], 
the area of origin of the six-vertex model with DWBC itself. 

The general procedure which we shall follow to build Fredholm determinant 
formula for the partition function has been suggested in [20]; in contrast to that 
paper, we shall however apply this procedure directly to the Hankel determinant 
representation (jl.4|) . In this way, the Hankel structure is preserved, and the 
integrability of the integral operator in the Fredholm determinant is ensured by 
construction. We therefore propose a simple factorization for the determinant 
of the Hankel matrix H appearing in (jl.4j) . This factorization, quite natural 
in the construction of a Fredholm determinant representation for the partition 
function, allows moreover to identify a core term in the factorized form of (|1.4J1 . 
with all others factor being trivial, and disappearing with a mere redefinition of 
vertex weights. The core term turns out to correspond exactly to the partition 
function for the six- vertex model when its i?-matrix is specialized to the ^(s^)- 
invariant .R-matrix [21], i.e. when its vertex weights are chosen accordingly with 
its underlying quantum group symmetry. This core term, which we shall denote 
as Zn gives rise to the Fredholm determinant, thanks to standard techniques 
from the theory of Orthogonal Polynomials. 

Together with this Fredholm determinant representation, we readily get for 
Zjy an equivalent representation as the ordinary determinant of an N x N sym- 
metric matrix, whose entries can be explicitly evaluated. From this last repre- 
sentation we propose a "reconstruction" formula which expresses the partition 
function as a trace of some quantum operator, which, in the large N limit, turns 
into the exponential of the boost operator for free fermions on a lattice. This 
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recalls analogous formulae which appears in the framework of corner transfer ma- 
trix and vertex operator approaches to integrable lattice models [8,22]. For iV 
finite, a corresponding "reconstruction" formula for Zn can also be written, as the 
trace (over the Fock space of N canonical fermions) of a product of exponentials 
of local operators. 



2 The factorization of the Hankel determinant 

In our way to build a Fredholm representation for the partition function Zn 
we shall follow the procedure suggested in [20], but applying it directly to the 
Hankel determinant representation (jl.4|) . As a consequence, the Hankel struc- 
ture is preserved, and the integrability of the integral operator in the Fredholm 
determinant is ensured by construction. In this section we shall discuss a specific 
factorization (a somewhat trivial one, but nevertheless basic for what follows) for 
the determinant of the Hankel matrix H appearing in (jl.4j) . 

All the different equivalent representations for the partition function derived 
in this paper stem essentially from the following simple observation: using the 
identity 

cot(A - 77) - cot (A + 77) (2.1) 



sin(A — 77) sin(A + 77) 
the matrix H can be naturally written as a difference of two matrices 

H = A_-A + , A ± = A\^ ±1 ± = A±77, (2.2) 

where the matrix A can be chosen to be 

A jfc = — [cot0-i] . (2.3) 

Our choice of the additive constant — i will be explained below; its role is to 
ensure invertibility of matrix ()2.3|) for any complex value of (p. The structure of 
(12. 2 j) suggests to factorize the determinant of H, for instance, as follows 

det N H = det N (A_) det N (I - AZ 1 A + ) . (2.4) 

The most evident consequence of such factorization is of course the relatively 
natural and straightforward emergence of a Fredholm determinant representation 
for the partition function, the Fredholm determinant being related to the second 
factor in (|2.4|) . The proposed factorization however suggests more: as we shall 
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discuss in detail in the following (see Sect. 4), a "reconstruction" formula can 
be deduced, which allows to represent the partition function as the trace of a 
suitably chosen quantum operator, in the spirit of corner transfer matrix and 
vertex operator approaches to integrable spin models. 

The analysis of factorization (|2.4|) relies essentially on the properties of matrix 
A and in particular of its determinant. To this purpose standard techniques 
relating Hankel matrices with orthogonal polynomials [23] will be exploited. Let 
us assume that entries of some generic N x N Hankel matrix A are given as 



A jk := / x j+k u(x) dx . (2.5) 



Let us moreover suppose that there exist a (complete) set of polynomials p n (x) 
orthonormal with respect to the measure u(x) dx : 



Pj(x)p k (x) u{x) dx = 5 jk . (2.6) 

Then, calling K n the leading coefficient ofp n (x), 

p n (x) = K n x n + ... , K n ^ , (2.7) 
the determinant of matrix (|2.5jl is simply given as 

JV-l 

det N A = Y[ Kj 2 . (2.8) 

3=0 

Of course, this formula turns out useful provided that the set of orthogonal poly- 
nomials associated to the measure u(x) dx can be identified. In the case of 
matrix H given by ()1.5|) appropriate polynomials are not available and the previ- 
ous scheme cannot be fulfilled for generic values of vertex weights. However, the 
matrix defined in equation ()2.3j) is much simpler and the scheme can be fulfilled 
explicitly. 

The entries of matrix ()2.3|) being periodic in Re 0, we may restrict ourselves 
to consider values of <fr varying over the vertical strip ^ Re0 < tt (with the 
point = excluded). In this region we may use 

1 — ™ dx ( 2 - 9 ) 

-oo 1 - e 7 ™ 

to write the entries of matrix A, Eq. ()2.3j) . in the form (|2.5j) with 

4>x 

u)( x ) = . (2.10) 

v ; 1 - e ra + i0 y J 
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The polynomials p n (x) (depending on as a parameter), associated to the matrix 
A, should therefore satisfy 



e 0x 



p j (x)p k (x) - _ ^ - - dx = 5 jk . (2.11) 

To identify the explicit form of these polynomials we shall now reexpress the 
orthogonality condition in such a way that the integration contour, though still 
being the real axis, has no singularity in its vicinity. This can be achieved by 
shifting the integration contour M — ^ IR — i and simultaneously relabelling the 
integration variable: x — > x — i. The orthogonality condition now reads 

/OO (f)X 
Pj (x - i) Pk (x - ^Y^dx = 8 jk . (2.12) 

Rewriting the weight function as 

1 rfi^Vf— (2.13) 



1 + e™ 2tt 

and comparing (|2.12|) with the orthogonality condition of Meixner-Pollaczek poly- 
nomials [24] 



1 r°° 

- / pW(x; 0) pW(x; <P) T(A - ix) r(A + ix) dx 
F(n + 2A) 



n!(2 sin^) 2A 

where 



5nm (2.14) 



PW(x;0) = ^^ 1 (^^ +il |l-e-*) 



(2.15) 



we readily identify the polynomials in question with nothing but the Meixner- 
Pollaczek polynomials, where the parameter A (not to be confused with the vari- 
able A entering the parametrization of the vertex weights) must be specialized to 
the value A = 1/2. 

Therefore, the polynomials p n (x) satisfying (|2.11|) are 

x + i 



p n (x) = e^ 2 P^ I <j>\ (2.16) 

or, due to (|2.15jl . more explicitly 

/ — n it /9 _. \ 

(2.17) 
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From this expression the highest coefficient K n of the polynomial p n {x) can be 
explicitly evaluated 



/2 (sin rh\ n + 1 / 2 



«n(^) = —f 1 > ( 2 - 18 ) 



nl 



and, due to formula (|2.8|) . we readily get for the determinant of the matrix A : 

i N<j) N-l 

det ^=7^V^ri(- ! ) 2 . (2.19) 



This last expression shows that if matrix A exists (that is if all its entries are 
finite) then it is invertible, since its determinant never vanishes (for finite 0). 
Factorization ()2.4|) therefore leads to an equivalent representation for the partition 
function, valid for all (non vanishing) values of </>±. 

Let us note that our choice of the constant (— i) in the definition of matrix 
(12. H)) can be easily explained now. Indeed, by considering a combination of (j2.19j) 
and its formal complex conjugate (i — > — i) one finds 



det 



dj+k ' jN 1 ™W) + "MW)^( nl) * m (22Q) 



■ (cot + a) 



The expression cos(iV0) + a sin(iV0) possesses zeroes in the complex plane of the 
variable unless a = ±i. Thus, by choosing a = — i (or, equivalently, a = i) the 
invertibility of matrix A is ensured. 

With (j2.19j) taken into account the partition function now reads 

Z N = [an e~ m - det N (I - AZ 1 A + ) . (2.21) 

It is worth emphasizing that the partition function is, in fact, described only by 
the last factor, the first two factors having a trivial meaning, and disappearing 
with a mere redefinition of the vertex weights. Indeed, the first factor in (J2.21j) 
can be seen as the result of a common prefactor in all weights, and we get rid 
of it by changing the overall normalization of the weights. The second one is a 
"boundary" factor specific of the DWBC choice for the six-vertex model; it can 
be removed by introducing a suitable asymmetry in weights W5 and W6, since any 
configuration contributing to the partition function (that is, satisfying both the 
"ice rule" and the DWBC), is such that the numbers of vertices of these two 
types satisfy the condition #w 6 — #w 5 = N, Hence, by choosing the weights to 
be 

wi = W2 = 1 , W3 = W4 = b/a , 

w 5 = e^^c/a , w 6 = e^^c/a , (2.22) 
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with a, b, and c given by ()1.3|) . the partition function reduces to the sole deter- 
minant in ()2.4|) . 

The choice of the weights ()2.22|) has a simple interpretation in terms of the 
six-vertex model i?-matrix, the matrix of local vertex states. In the standard 
notations the choice of the weights in the form ()2.22|) corresponds to the following 
i?-matrix 



/l 



R{y) 



0{u) 
e" iy 7(z/) 



e iiy 7 (z/) 

Ply) 



\ 



V 



(2.23) 



where 



Ply) 



Sill V 



sin(t/ + 2rj) 



7(1/) 



sin2r/ 



sin(v + 2rj) 



v 



A 



(2.24) 



In this expression one easily recognizes the [^(s^-invariant i?-matrix [21] which 
is, moreover, normalized to satisfy the "unitarity condition" : 



R(u) (PR(-v)P) 



(2.25) 



with the permutation operator P := R(0) . Thus, the factorization (J2.4j) of the 
initial Hankel determinant leads naturally to the representation of the partition 
function just as the sole determinant in (|2.2ip . provided that the vertex weights 
are chosen according to the underlying quantum group symmetry of the model. 
From now on we shall assume this choice, equation (|2.22jl . for the vertex weights, 
denoting the corresponding partition function as Z N , given as 



det N (I- AZ A 



(2.26) 



It is worth mentioning here that almost everything above can also be extended 
to the case of the determinant formula of Ref. [2] for the partition function of the 
inhomogeneous model. Without giving any detail, let us just emphasize that pre- 
vious considerations become even more transparent, since in that (more general) 
case matrix A is a Cauchy matrix. However, the homogeneous model possess 
the further interesting property that the partition function can be expressed as 
the Fredholm determinant of an integral operator of integrable type, in the sense 
paper of [18]. This will be the subject of the next Section. 
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3 The partition function as a Fredholm determinant 



We shall now focus our attention on the determinant in ([2.26)1 . First of all we 
need to build the entries of matrix Az 1 = A~ 1 \^ )=( j ) _. This task can be achieved 
straightforwardly by borrowing standard techniques from the theory of orthog- 
onal polynomials [20,23]. Let us recall that, indeed, once the set of orthogonal 
polynomials associated to a given Hankel matrix A with entries (|2.5jl is known, 
the entries of the corresponding inverse matrix A~ l can be evaluated from the 
function 

N-l 

IC N (x,y) = ^2p n (x)p n (y) . (3.1) 

n=0 

simply in terms of partial derivatives: 

The proof is based on the fact that function Kn{x, y) is the kernel of an integral 
operator, with respect to the measure u>(x)dx . This operator, by construction, 
projects over the subspace of polynomials of order less than N, and therefore acts 



as identity operator on monomials l,x, 
JC N (x,y)y m Lj(x) dx 



x 



,x 



m 



n-i . 



0,1, 



N-l 



(3.3) 



With ()3.3|) taken into account, equation ()3.2|) can be verified directly: 



N-l 



N-l 



1 1 d k 



k=0 



k=0 

1 d j 



z k+m u(z) dz 



j\ dxi 
1 d j 



j\ dxi k\ dy 

lC N (x, z) z m uj(z) dz 



x=0, y=0 



x=0 



a;" 



2=0 



j\ dxi 

= S jm . (3.4) 

In the case of matrix A ([2.3)1 the entries of the inverse matrix A" 1 are given by 
(13.2)1 with the function K,N{x,y) build from formula ()3.1[) in terms of orthogonal 
polynomials p n (x), defined in Eq. (J2.17)) . 

Generalizing the simple computation above one can evaluate the entries of 
the product matrix AZ 1 A + : 



J"' 



1 d j 
j\ dxi 



JC^(x,z) z m u + (z) dz 



(3.5) 



x=0 
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Here the + or — superscripts denote the dependence on the variables <j>± respec- 
tively. The determinant in ()2.26|) can now be transformed as follows 

In detjv(I - AZ 1 A + ) 

oo 



n=l 
oo 



— 1 f°° f°° 

n=l 71 1=1 

oo 



n=l 



(3.6) 



(3.7) 



(3.8) 



= lndet(T-V7v) 
where Vat is the integral operator on the real axis with kernel 

V N ( x ,y) = Kirfay) uj + (y). 
Using the Christoffel-Darboux identity, this kernel can be written as 

V N {x,y) = — iu^{y) 

k n x y 

rendering the integrability (in the sense of Ref. [18]) of integral operator Vat 
manifest. Integral kernels of the form ()3.8|) are also known under the name of 
correlation kernels since they arise in expressions for eigenvalues correlation func- 
tions in the theory of random matrices. For discussion of the role and importance 
of this special class of integral operators in connection with correlation functions 
of integrable models and with the theory of random matrix, see [19,25]. 

Taking into account formulae ()2.10|) . (|2.16|) . (|2.17|) and (J2.18|) we therefore 
get for the partition function of the model, with the choice ()2.22|) for its vertex 
weights, the following Fredholm determinant representation: 



Z N = det(J- Vat) 
where the kernel may be most explicitely written as 



(3.9) 



Vn(x,v) 



2 -Pi 
— 2F1 



-N, ix/2 
1 

-N + l, ix/2 
1 

jy e 2iN</>_ e <j} + y 



l_ e -^-j 2jPl 
1 - e" 2 ^ 



N + 1, iy/2 
1 

-N, iy/2 
1 



1 - e" 2 ^ 
l_e- 2i ^- 



x 



x 



y 1 - e*y + iO 



(3.10) 
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This representation is valid for < Re0 + < tt and arbitrary complex </>_. Let us 
moreover underline that, contrarily to the original Hankel determinant formula, 
the parameter N can be extended here from the set of positive integers to the 
whole complex plane. 

We shall now discuss alternative forms for the representation we have just 
obtained. First, it should be noted that by shifting the integration contour and 
relabelling the integration variables in each integral in ()3.6|) this result can also 
be put in the form 

Z N = det (J - (W N ) , C = e^— , (3.11) 

with (to shorten the formulae we shall use more compact notations in terms of 
polynomials) 

W N {x, y) = N 

x-y 

X 1 + e 2 ^ ' (y3 ' 12 ^ 

We see here that as a matter of fact the kernel is a real-valued function for real 
(j)± (since Meixner-Pollaczek polynomials P„ (x; <fi) are real- valued for real x and 
0). 

The parametrization of the weights in the form (jl.3j) with real A and rj (hence, 
real </>±) is typical for the so-called disordered phase of the model [8]. Thus, the 
just obtained representation for the partition function, even if valid for arbitrary 
choice of vertex weights, can be regarded as "adapted" to the disordered phase. 
The other two physical regimes are ferroelectric and antiferroelectric and they 
can be obtained by choosing <fi± to be purely imaginary. 

To obtain corresponding representations it is sufficient to note that in the case 
of complex + the integrals in (|3.6J) can be evaluated by closing the integration 
contours upwards (downwards) in the complex plane of the variables x 1; . . . , x n if 
Im + > (if Im + < 0). As a result, each integral is given by a sum of residues at 
simple poles of the function u{x) lying in the upper (respectively lower) complex 
half-plane. Considering the case of Im0 + > (the case of Im0 < leading to 
an essentially equivalent result) we obtain the following representation for the 
partition function: 



■IN 



det(J-Vjv) (3.13) 
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where Vn is the integral operator with discrete kernel 
V N (x, y) = -^M N (x; 1, e~^-)M N _ 1 (y; 1, e" 2 *-) 

-M N ^{x;l,e-^-)M N (y;l,e- 2 f-)\— e" 2 ^ 

(3.14) 

whose "integration" variables x, y take nonnegative integer values; for x = y the 
kernel is to be understood in the sense of the Cristoffel-Darboux identity; in the 
last formula, the standard notation for Meixner polynomials 



M n (x;(3,c) = 2 F 1 



-n, —x 



1 - - ) (3.15) 



c 



has been used [24]. The "tilded" variables are defined as ± = i</>± . The 
representation ()3.13|) is valid for Re <p + > and arbitrary complex </>_ . The real 
values of <j>± correspond to ferroelectric (0± > 0) and antiferroelectric (0 + > 0, 
4>- < 0) phases of the model. 

It is to be mentioned that representation (|3.11j) can be also obtained directly 
by employing the formula 

oo 

coth0= 2^e" 2 ^ , Re0>O (3.16) 

x=0 

instead of (|2.9|) and repeating all considerations of the previous Section. The 
appearance of Meixner polynomials is then quite obvious since for |c| < 1 they 
are subject to the orthogonality condition 

V Mfa (3, c)M k (x; (3, c)^% c x = 5 jk (3.17) 
^— i x\ 1 — c 

and upon setting (3 = 1 and c = e -2 * the identification of the proper set of 
orthogonal polynomials is achieved. 

We conclude this Section by considering the so-called rational parametrization 
of the vertex weights, which corresponds to the case in which a, b and c are 
restricted by the condition a ± b = c . This regime can be obtained through a 
suitable limit from vertex weights (jl.3j) . Namely, depending on the choice of the 
sign in this restriction, one should just substitute A, rj —> eX, erj (for plus sign) or 
by A, r] — > 7r/2 — eA, 7r/2 — er] (for minus sign) in (|1.3|) and take the limit e — > 0, 
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after renormalization of the weights by a factor 1/e. We should of course recover 
in this limit the corresponding result of Ref . [20] . 

The Hankel determinant formula for the partition function of the model with 
rational weights is simply given by formula (|1-4|) with the rational functions A, rj 
instead of sine functions: sin(*) — > (*). The Fredholm determinant representation 
for the partition function in this case involves Laguerre polynomials; the "rational 
limit" of the kernel ()3.12j) can be easily found using 

hmP^x/e;^) = L„,(-20x), ^oa—-— = e^-y), (3.18) 

where L n (x) is Laguerre polynomial, and 9(x) is Heaviside step-function. The 
same result can be obtained from ()3.14j) . too, with the discrete measure turning 
into the continuous one in the standard way, when performing the rational limit. 
Explicitly, the partition function in the rational case is given by the Fredholm 
determinant of the integral operator on the real positive half-axis, defined by the 
kernel 

V„(x, y) = _ N ^x)L N ^y) - L N ^x)L N m ^ 

x-y 

where £ = The result of Ref. [20] is thus reproduced in the rational limit, 

which corresponds there to the case 5 = 1. 



4 The finite-size determinant representation 

The determinant in ()2.26|) can be also written as that of some N x N symmetric 
matrix whose entries are simply connected with the kernel of the integral operator 
in the Fredholm determinant representation. One can insert (|3.1j) in the second 
line of and obtain 



lndetjv(7 - AZ 1 A + ) 

= ~y^~ / • • • / ^n(x 1 ,x 2 ))C^(x2, x 3 )... IC N {x n , xi) JJ uj + (xi)dxi 

1 " J —OO J —OO i i 



n=l 



-oo 

oo JV-1 



£ ^ _ k 2 k 3 ■ ■ ■ v k n k 1 

n=l k\,...,k n =0 



= lndet N (I -V) (4.1) 
where entries of the matrix V are 

POO 

v jk= p~(x)pZ(x)uj + (x)dx . (4.2) 
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Clearly, this identity is a consequence of the fact that the matrices AZ 1 A + and 
V are related by some similarity transformation. 

Hence, the partition function admits also the representation 

Z N = detjv(/- V). (4.3) 

It turns out that in contrast to entries of matrix H, entering the original Hankel 
determinant representation (|1.4|) . those of the matrix V can be computed ex- 
plicitly. It is convenient (in analogy with (|3.11Jl ) to introduce the matrix W by 



V = (W , C = e^--*+> = e" 2i7? . (4.4) 
Entries of W are defined by the formula 



„2:r© . 

H/, fc = 2sin0_ / pf/ 2 \ x -^)P^ 2 \ x - (t) ^ Y -^-dx. (4.5) 



As shown in detail in the Appendix, the integral can be evaluated in closed form. 
Introducing, for the sake of convenience, the notations 

sin0_ sin(A — rj) sin(0 + — </>_) sin 2^ 

p = ~^~ir = ■ ( \ ; — v ' 7 = — ~a — = • T\ ; — V ' 

sm0 + sm(A + 7/J sin0+ sm( A + r\) 

(note that (3 and 7 are precisely the quantities entering the £/q (5 [^-invariant 
.R-matrix, see eqs. ()2.23|) and (|2.24|) . (3 = /3(is), 7 = 7(1^)) the result reads 



"•—-§(;)(:)(?) 



n=0 

?2 



-j, ~k 
1 



EL 

sy2 



(4.7) 



It is worth to remind that the representation ()4.3)1 is for the partition function 
with the weights ()2.22j) : the original normalization can be achieved by multiplying 
the RHS of (JOD by the factor [sin(A + rj)] N2 e^ x ~^ N . 

The second expression in (I4.7J1 shows that the entries of matrix W are, in fact, 
Meixner polynomials, see (j3.15j) . On the other hand, from the first expression in 
(14.7)1 follows that matrix W can be written as a product of much simpler matrices. 
Indeed, introducing N x N matrices 

(J+)nm ^ $n— l,m 

(Jo) nm (n + 1/2) 5 n 

,711 

{J-)nm = (n + 1) <Wl, m ( 4 -8) 



14 



and taking into account that 

( ex P{7^ + })„ m = 7"- m (^), J- = (J + f, (4-9) 
one obtains 

W = exp{ 7 J + } (3 2J ° exp{ 7 J_} . (4.10) 

Thus, the matrix W is defined by its Gauss decomposition, i.e., as a product of 
a lower-triangular, a diagonal, and an upper-triangular matrices, respectively. 

Matrix W simplifies considerably in the limit iV — > oo. This is a consequence 
of the fact that the semi-infinite dimensional matrices J±£ with entries (|4.8jl 
satisfy su(l, 1) algebra commutation relations 

[J_,J + ] = 2J , [J±,J ] = TJ±- (4.H) 

Standard techniques (widely used for instance in the theory of generalized coher- 
ent states [26]) can now be applied to (|4.10j) . with the following result 

W = exp{2r]K} (N = oo) , (4.12) 

where the semi-infinite dimensional matrix K is 

K = — ( J_ + J + - 2 cos v J ) . (4.13) 
sm v v 

Here (and below) we use "mixed" set of variables, v = 0_ = A — r\ and rj (indeed 
the most natural ones for the U q (s\>2) -invariant i?-matrix, see Eqs. ()2.23|) - (j2.25)) ). 

The matrix K is the Jacobi matrix corresponding to the recurrence rela- 
tion of Meixner-Pollaczek, Laguerre or Meixner polynomials depending whether 
| cosz/| < 1, | cosz/| = 1 or | cosz^| > 1, respectively. These polynomials are there- 
fore eigenf unctions of K with eigenvalue x, and the spectrum of K is thus given 
by the support of the measure appearing in the orthogonality condition of the 
corresponding polynomials, which is R, M + or Z + , respectively (i.e., it is continu- 
ous in the first and second case and discrete in the third one). From this point of 
view the presence of these same polynomials in the kernel of the integral operator 
in the Fredholm determinant representations is obvious: the integral operator is 
just the semi-infinite dimensional matrix 11^ exp{2rjK } re-expressed in the basis 
of the eigenfunctions of K\ here the iV-dimensional projector II at is the matrix 
with its first iV diagonal entries equal to one, and all others entries equal to zero. 

Another point which is to be discussed is that formula (j4.12j) tells us that in 
case N = 00 , matrix W (and hence V) is just the exponential of some "simple" 
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matrix. The analogy which come into mind after looking at formulae ()4.3|) and 
(j4.12j) is with a typical result of calculation of traces of certain class of operators 
acting over the Fock space of N fermions. Namely, given a matrix A, and the 
corresponding quantum operator A, bilinear in canonical fermion operators, built 
from matrix A as follows 



JV-l 



A — ^ ] cl l+1 A nm c m+1 , c^c m + c n $ m — 5 nm (4-14) 

n,m=0 

it is well-known [27] that 

Tr[expA] = det7v(/ + expA) . (4.15) 

To make a connection with formula ()4.3|) . it is to be mentioned that the minus sign 
in RHS in (|4.15jl can be acquired by considering supertraces instead of traces; the 
supertrace is defined as Str[*] := Tr [(— 1)-^ * ] , where N denotes the fermion 
number operator: N = Yl n =i ^n^n • However, since the fermion number operator 
commutes with any operator of the form ([4. 14)1 . one can consider just traces when 
dealing with the exponentials of such operators. Taking into account this we can 
therefore write the following trace formula for the partition function 

Zn = Tr exp {2vLn}\ (4.16) 

with 

L N = K N + i^^N , (4.17) 

where the operator Kn is built from the matrix 

K N = — In W . (4.18) 

Note that in contrast to the entries of matrix W those of matrix (and hence 
of Ljv) essentially depend on iV as indicated by the notations. Operator K^, 
though bilinear in fermions, is essentially non-local, thus making the entries of 
matrix for finite N rather complicated quantities. However, formula (|4.12j) 
implies that matrix simplifies considerably in the case N = oo, and the 
corresponding operator K := becomes a local one. Explicitly we obtain 

oo 

K = J2 n] ^n,n+l, (4.19) 

n=l 
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where 



// 



n,n+l 



sm v 



St 



cosz/(c+c„ + c+ +1 c n+1 ) 



(4.20) 



We recognize in this expression for operator .fiT what is known in the literature 
as a boost, or ladder, operator. To be precise K is the positive half- axis part 
K = flW of the total boost operator B = B^+B^ + \ with [B( + ),.BH] = 0, see 
[28]. Operator B is the boost operator for the model described by the following 
Hamiltonian 



if 



1 



sm i/ 



// - 2 cot uN 



(4.21) 



where Hq is the hopping term. From this formula it is clear that u, the spectral 
parameter in the six-vertex model, is here playing the role of a chemical potential 
rather than a coupling constant. Nevertheless, all the construction above has 
led us to a description which is quite analogous to what we have in the corner 
transfer matrix formalism [8], or to a more general extent, in the vertex operator 
approach to integrable models [22] and angular quantization method in integrable 
quantum field theory (see, e.g., [29] an references therein). 

As previously explained, in the case of finite iV operator is nonlocal 
(though free-fermionic). One can however still write the partition function in 
terms of the product of exponentials of local operators. As a generalization of 
()4.15|) the following is also valid 



Tr 



Y\ exp Ai = det N ( I + exp A { j 



(4.22) 



where operators are constructed out of matrices Ai through formula (|4.14|) . We 
can therefore write the following "reconstruction" trace formula for the partition 
function 



with 



Tr 



JV-l 



e- 2iv ) N exp^J+j (3 ZJ ° exp{ 7 J_} 



2J n 



N 



J+ = ^2 nC n+l C n , Jo = J^t 71 ~ X / 2 ) C « C « ' 3 = & 



n=l 



n=l 



(4.23) 



(4.24) 



The reconstruction formula ()4.23j) expresses the partition function of the six- 
vertex model as the trace of a quantum operator built out of fermions, or, modulo 
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Jordan- Wigner transformation, of spin-1/2 operators. Note that, in fact, we have 
derived these expressions for the partition function starting from the Hankel de- 
terminant representation, previously obtained within the Quantum Inverse Scat- 
tering Method. The question which arises from the considerations presented here, 
is how representation ()4.23|) in terms of fermions (or spins), which are to be in- 
terpreted as effective degrees of freedom, could be extracted directly from the 
basic definition of the model in terms of vertex configurations. An answer to this 
question might suggest alternative approaches to the open problem of calculation 
of correlation functions for the model. 
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Appendix 



P^\x- t) P«(i; u) |r(A + ix)| 2 e (2< ^ dx (A.l) 



The integral defining the entries of matrix W in (|4.5jl is a particular case of the 
integral 

j£>( Tjtl ;;0 = J- 
271 

where <fi G (0, it) and A is assumed to be real and positive, A > 0. These restriction 
are important for convergence of the integral. Our aim here is to prove that for 
arbitrary r and u the quantity Imn(r, u; 4>) has the following expression 



r(2A + n)r(2A + m) 
(2sin0) 2A r(2A)n!m! 
— n, —m 
2A 



sm r 



x 2 F X 



sm (p 
sm r sin u 



sin w 



sin i 



sin(r — (p) sin(co> — 



(A.2) 
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Equality (jA.2|) is a consequence of the orthogonality condition for the Meixner- 
Pollaczek polynomials and of the identity 

i n—k / . \k 



p(A)/ \ _ V- T{n + 2\) [sin(r-0)]" (sinr) (A) 

(A.3) 



, Q T(k + 2X)(n-k)\ (sin 



which can be viewed as the extension of an analogous formula for Laguerre poly- 
nomials (see, e.g., Ref. [30], §10.12, Eqn. (40)). Indeed, using identity (jA.3|) 
for both polynomials under the integration sign and employing the orthogonality 
condition (J2.14|) one immediately obtains 

(A ), , [sin(r-0)] n [sin(^- 



2 2A (sin0)™+ m + 2A 

min(n,m) 

X 



^ r(n + 2A)T(m + 2A) 



T(k + 2X)(n~k)\(m-k)\k\ 



sin r sin u 



sin r — 0) sin w — 



(A.4) 

Rewriting the finite sum here as a truncated hypergeometric series results in 
expression (|A.2J) . 

Let us prove now the identity (|A.3|) . To simplify as much as possible the 
combinatorics, let us consider the three-term relation satisfied by the Meixner- 
Pollaczek polynomials: 

(n + I) P^x; <p) - [2x sin + 2(n + A) cos^] P^ A) (x; 0) 

+ (n + 2A-l)i^ A) 1 (x;0) = 0. (A.5) 

Let us define 

Snix-A) = (sin0)~ n pW(x;0) . (A.6) 

Recurrence relation (|A.5|) takes the form 

(n + 1) S n+ i(x; 4>) - 2(n + A) cot</> S*„(x; 0) 

+ (n + 2A- l)[l + (cot0) 2 ] S„_i(x;0) = 2a; S„(:r; 0) . (A.7) 

From this relation it is clear that S n (x; <fi) depends on <p only through cot <f), and, 
moreover, it is a polynomial of order n in cot <fi. It is useful to consider the case 
(p = 7r/2 so that cot = 0. Denoting S n (x) = S n (x; ir/2), in this case one has 

in + 1) S n+1 (x) + (n + 2\- l)S n ^(x) = 2xS n {x). (A.8) 
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Let us consider the semi-infinite dimensional matrices 



( J °)nm = ( n + X ) 5 n,m 

{J + ) nm = (n + 2X-l)5 n ^ m (A.9) 

where, as in the main text of the paper, n, m 6 {0,1,2,...}. These matrices 
obey su(l, 1) algebra commutation relations (|4.11|) . Modulo a diagonal similarity 
transformation, they correspond to the standard matrix realization of the positive 
discrete representation X>W(A) of su(l, 1). In terms of matrices (|A.9J) relations 
fTOjl and (fA~8j) read 

[ J- - 2 cot J + (l + (cot 0) 2 ) J+] £(z; 0) = 2x 0) 

[J_ + J+] £(a?) = 2x £(z) , (A. 10) 

where n-th component of S(x,<f)) is just S n (x,<f)), and S'(x) = S(x,tt/2). Using 
the commutation relations of su(l, 1) algebra it can be straightforwardly checked 
that the following relation is valid 

exp{aJ + }[j_ + J + ] = [j_ - 2aJ + (l + a 2 ) J + ] exp{aJ + }. (A.ll) 
Thus, by identifying a = cot we obtain 

S(x;cj)) = exp{cot0 J+} S(x) . (A. 12) 

This is a key identity to prove relation ()A.3j) . Indeed relation (|A.12J) implies that 



— * — * 

S(x; t) = exp{(cot r — cot 0) J+} S(x; 0) . (A. 13) 

Then, taking into account that 

(exp{aJ + }) = - r (rc + 2A ) a " , A -m 

one obtains 

S n {x-r) = X) r? ^9^ 2A) M (cotr-cot0r- m g m (g;0). (A.15) 
^ o r(m + 2A) (n-m)! 

Rewriting the last equation in terms of the standard Meixner-Pollaczek polyno- 
mials, see flA.6|) . one arrives finally to identity ()A.3|) . which is thus proved. 
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